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Abstract 

Thermodynamic properties of ideal Bose gas trapped in an external generic power law potential are investigated 
systematically from the grand thermodynamic potential in d dimensional space. The most general conditions for Bose- 
Einstein condensate and the discontinuous conditions of heat capacity at the critical temperature in presence of generic 
power law potential are presented in this manuscript. The dependence of the physical quantities on external potential, 
particle characteristics and space dimensionality are discussed. The more general results obtained in this paper 
presents a unified illustration of Bose-Einstein condensation of ideal Bose systems as they reduces to the expressions 
and conclusions available in the literature with appropiate choice of power law exponent. 

Introduction 


Many attempts have been made to understand the phenomenon of Bose-Einstein condensation (BEC), after the demon- 
stration by Einstein that there is a possibility of condensation of free bosons mm- Since then the bulk behavior of bose 
gas are investigated by many authors [SI jHISl EllTllHl IB HB [H [IB • An increasing attraction towards this subject is ob- 
C served, after the achievement to create BEC in magnetically trapped alkali gases [Ml[13[IS]. As a result, different studies 
0 analyzing effects of external potentials [TT) [IB [IB HB HD: relativistic effect [TTl HB [S], space dimensionality [T71 HB UBi 
i__icriticality near the transition point [SB [21], casimir effect [SB [13 [2S] of the bosonic system have been performed. Moreover, 
the recent experimental observation [SB of condensation of massless bosons makes the research area even more challenging. 


In a real system, of course interaction between particles do exist. But taking it into account makes the problem difficult to 
solve analytically. Interactions between bosons and their effect on critical temperature Tc and fraction of condensation ^ 
are understood for experimental situation, which are few percent or less when the density of gas is low [IB. Thus it is well 
r approximated that the Bose gas of low density can be treated as in ideal bose gas. Neveretheless, to understand the effect 
of interactions and to retain the essential physics, we approximately represent the real system by non interacting particles 
(N in the presence of an external potential such as harmonic potential [IB, toroidal potential[TB, power law potential [171l22] . 
mean field potential[23. The constrained role of external potential for atomic gases do change the performence of gases. 
Thus trapped atomic gases provide the opportunity to manipluate the quantum statistical effects. An excellent study by 
• • Luca Salanasich [17jon Bose gas confined in a trap, described by a symmetric (isotropic) power law potential {U = Ar"^) 
. ^ d-dimensional space presents a very important relation whether the system will exhibit condensation or not. Salanasich 
draw a conclusion from the argument of number of particles that, BEC exists if f -k ^ > 1, where n is the exponent 
of symmetric power law potential. Also, from this relation one can say BEC is possible in trapped Bose gas under this 
^ type of potential in d = 2 with appropiate choice of n, which was not the case for ideal free Bose gas. So, it would be 
interesting to investigate the relation when the potential is not symmetric in general. A point to note, BEC for free 
bosons at d = 3 is a first order phase transition[3 El [B- But the order of phase transition for the trapped bosons in 
arbitrary is yet to be examined. 


Lot of efforts have been made in order to understand dimensional dependence of phase transition [^ l9l ITTl IT^ [TTl [2(711^ 1^ . 
Thermodynamic as well as other properties are calculated in d-dimension for bosons as well as Fermions [301 UTl 132] and 
interesting conclusions such as equivalence relation between the Bose and Fermi gases in two dimensions is established|33j. 
In this report we intend to calculate all the thermodynamic properties of ideal Bose system in d dimension, under generic 
power law potential (not essentially symmetric). In order to do this, we have first calculated the grand thermodynamic 
potential. And then from grand potential, we have calculated the thermodynamic quantities such as internal energy E, 
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entropy S, pressure P, number of particle N, helmholtz free energy A, isothermal compressibilty /c, specific heat at con¬ 
stant volume Cv and pressure Cp, their ratio 7 = From the thermodynamic quantities we have derived important 
results regarding the Bose system such as condensate fraction, specific heat jump, latent heat of condensation, critical 
temperature, equation of state as well as general criteria for existance of BEC for trapped boson under generic power law 
potential. Beside this, from the Clapaeyron equation it has been shown that, like the free ideal bosons in d=3[7ll9], BEC 

for trapped bosons is also a first order phase transition in any dimension. A point to note that, in the hamiltonian instead 
2 

of ^ type kinetic part, we have took ap^, where p is momentum and a is constant. And thus making the important 
conclusions in more generalized way and obviously the results obtained here provide a unified description of behavior of 
trapped Bose gases, from which many results in the literature can be derived. 

The report is organized in the following way. The density of states and grand potential are calculated in section 2. 
Section 3 is devoted to investigate the thermodynamic quatities and to check the order of phase transition. Results and 
discussions are presented in section 4. The report is concluded in section 5. 


2 Density of States and grand potential of bose gas under generic power 
law potential in d dimension 

Let us consider the ideal Bose gas in a confining external potential in a d-dimensional space with energy spectrum. 


d 

e(p,a;i) = V+^Ci| —I"-’ (1) 

Where, b, I, ai, Ci^ rii are all postive constants, p is the momentum and Xi is the i th component of coordinate of a 
particle. Here, Ci, Ui and rii determines the depth and confinement power of the potential. With I = 2, b = ^ one can 
get the energy spectrum of the hamiltonian used in the literatures [amiin]. For the free system all rii —> 00 . 


Density of states can be obtained from the following formula, 


P(e) 



e{p,r)) 


Now, from Eq. (1), the density of states. 


p(e) = B 
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(3) 


where. 


gVdCd A r(;^ + 1) 


(4) 


Here, Cd = rp^f+T)’ ^ degeneracy factor, Vd = volume of an d-dimensional rectangular 

whose i-th side has length 2ai and r(l) = dxx^~^e~^ is the gamma function. 


The grand potential for the Bose system. 


q = — ln{l — zexp{—j3e)) 


(5) 


/3 = ^, where k being the Boltzmann Constant and z = exp(/3/r) is the fugacity, where fi being the chemical potential. 
In experiments with trapped gases, thermal energies far exceed the level soacingflb]. So, using the Thomas-Fermi 
semiclassical approximation [53] and re-writing the previous equation. 


q = qo- 


p{e)ln{l — zexp[—j3e)) 


( 6 ) 
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( 7 ) 


So, using the density of states of Eq. (3) we finally get the grand potential, 


q = qo + Sr(y + l){kTYg^+i{z) 

where, q^ = —/n(l — X = j + Ylt=i ;r 9i{^) is the Bose function. Where, 

/•oo ^i-i °° 

9i{z)= - T='I2~I 

Jo z -1 ^ 

Point to note when / > 1, as z —>■ 1 (in condensed phase), 

OO -| 

»(i) = i:? = cm 

r—1 

3 Thermodynamics of bose gas under generic power law potential in d 
dimension 

3.1 Number of Particle 

The number of particles N can be obtained. 


dq 

N = z{ — )py = No + 


gC^T{^ + l)VaY{Un^ + ^) 


Here, No = is the ground state occupation number. 

Now, defining, 

hT ^ 

= Pdn(—)'^”*r(- + i) 

rz 71 z 


(kT)^g^(z) 


(9) 


A' = 


Z =1 

hbr dl2 + 1. 


TTs (kT) ? d/l + 1 ■ 
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It is noteworthy. 


lim Vj = Vd 


lim X = 7 

7 i^OO t 


lim A' = An 
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( 12 ) 

(13) 

(14) 


So, if we choose I = 2 and b = -^ from Eq. (14) we get Aq = which is the thermal wavelength of nonrelativis- 

tic free massive boson. However, it should be noted that, when I ^ 2, X' then depends on dimension. With d = 3 and 
d = 2, thermal wavelength of photons are respectively { 2 -ny%kT can be obtained from from Eq. (11) 

by choosing b = c, where c being the velocity of light. So, from the definition of X' with more general energy spectrum, 
one can reproduce the thermal wavelength of of both massive and massless boson. 


The number of particle equation is then written as, 

N-No = g^gAz) ( 15 ) 

A 

With I = 2, a = all rii —> oo the number of particle equation for free massive boson in d dimesnsional space can be 
obtained from Eq. (12), (13), (14), (15), 

N-No=g^g.iz) (16) 

Ao 
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which is exactly like Ziff[5], which gives the exact equation for number of particles at d = spin. 


Now turning our concentration to find the critical temperature, we know T —> Tc, n —>■ 0 or z —> 1. So, from 
equation (9), the critical density is, 


Nr. = 


9CnT{j + i)uUn^ + i) 




llrii 


{kT)K{x) 


(17) 


And the critical temperature, 

^ 1 Ik 

’^9Cnnj +mutiny+ mx) 

Also from equation (17) and (18), we can obtain, the fraction of condensation, 


(18) 


N 


T 


X 


For free massive boson, the two above equations become like below as in Ref. [3 [25] 


(19) 


V^TTmk^ Cid/2) 

^ _ 1 _ 

N ~ ^tJ 


( 20 ) 

( 21 ) 


From Eq. (20), it is seen if d —> oo, = ( ^ 2 ^ which we can say as d approches infinity, the critical 

temperature picks a non-zero value instead of being zero m- Also Eq. (20) and (21) exactly reproduce the critical 
temperature and condensate fraction at d = 3[6l|7lli|9]. 


Now, from equation (15), one can find the most general criterion for BEC to take place. It is seen that whether 
BEC to take place or not depends on y. Eq. (15) dictates BEC to take place when, 




( 22 ) 


this is satisfied. This equation shows that the conditions for BEC to take place depend both on dimension d as well as 
on kinematics characteristic I and exponent of power law rii. When the potential is symmetric, i.e., ni = n 2 = ■■■ = rii = 
...Ud = n, Eq. (22) becomes, 


d 

I 



(23) 


Now, putting Z = 2 for massive boson in the above equation we achieve the criterion for BEC to take place in case of 
symmetric potential which is same as Ref. m 


d 

2 



In case of free massive boson, all Ui —> oo, so Eq. (24) becomes. 


(24) 


d 

2 


> 1 


(25) 


So, for free massive boson, we get the usual criterion that d should be greater than 2 for BEC to take place, consistent 
with the result obtained in Ref. [3135]. 


3.2 Internal Energy 

Erom the Grand Canonical Ensemble internal energy E is, 

F = 

gC„r(f+ i)Edntir(; 




1 ), 




( 26 ) 
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So, the internal energy below and greater than the critical temperature 


E = 


NkTx<^ 

{N - No)kTx^^ 


,T > Te 
,T<Te 


In case of free massive bosons {I = 2, b = ^ and rn —>■ oo), one can find Eq. (27) reduces to,. 


E = 


2 gd/2{z) 


,r > Tc 
,r< r. 


which is in accordance with ZifF[3], also produces the exact expression of if for d = 3[S1[7]. 


(27) 


(28) 


Now as T » Tc, from Eq. (27) it is seen, the internal energy becomes, E = NkTx- For free massive bosons it is 
E = ^NkT, which becomes ^NkT, when d = 3, thus E approaches the classical value at high temperature. 


3.3 Entropy 

The entropy S can be obtained from Grand Canonical Ensemble, 

^ = krA^y-Nklnz + kq 

al 

Again the entropy below and greater than the critical temperature 


5 = 


l)5x+i(^)-Inz] ,r>Tc 

(iV-iVo)fc^(x+l)C(x + l) ,T<Te 


(29) 


From the above it is easy to see S' = 0 as T —> 0, in accordance with the third law of thermodynamics. As before, for 
free massive bosons (/ = 2, 6 = ^ and rii —oo), with the help of Eq. (12)-(14) one can find Eq. (29) reduces to, 

( Nk[j^il + l)gd^,iz)-lnz] ,T > T, , , 

1 (iv-iVo)fcfJ(f+ i)C(i + i) ,T<T, ^ ’ 

which is in exactly like Ziff[9]. Again at d = 3 Eq. (30) reduces to same expression for entropy as Ref.[6l[7] 


3.4 Helmholtz Free Energy 

From the Grand Canonical Ensemble we get the expression of Helmholtz Free Energy, 

A = -kTq + NkTlnz (31) 


Now from the expression of grand potential we obtain the equation of Helmholtz Free Energy below and above 


A = 


NkT 

NkT 


9x+i(A 

SxiA 

C(x+1) 

C(x) 


,T>Tc 

,T<Te 


(32) 


In case of free massive boson the above expression reduces like below. 


A 

NkT 


+ 

(z) 
C(| + l) 

~cW 


+ \nz 
+ lnz 


,T>Tc 

,T<T, 


(33) 


exactly like Ziff[2. Now, for d = 3, the above equation produces the exact expression for Helmholtz free EnergyjBJ [7]. 
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3.5 Pressure 

Rewriting equation (9) stating the number of particles, 


N -N, 


0 


N-No g , , 




Now a very well known expression for the nonrelativistic d—dimensional ideal free Bose gas[S], 

N-No g 

= Jd9d/2{z) 


Vd 


Comparing the above equations, we can say is a more generalized extension of Vd- It represents the effect of external 
potential on the performence of trapped bosons. Calling the effective volume the grand potential can be rewritten as. 


q = 


gV^ 

qo + g^gx+i{z) 


(34) 


So, the effective pressure 


Also can be rewritten as. 


P' 


-( —) 


^gx+i(z) ,T>T, 

^C(x + i) ,t<t. 


^ I dq , 


NkT gx+i(z) 

Vd 3xV) 
{N-Nq) ,rr. ax+l) 
Cl Cix) 


,r > Te 

,T<T, 


(35) 


(36) 


The above equation is very general equation of state. It is the equation of state for any dimensionality d, any dispersion 
relation of the form (oc p'^) having any form of generic power law trap and obviously it is expected that it will reproduce 
the special case of free system. For free system the equation (36) becomes. 


P=h—) = 

rdVd’ 


NkT gd/2+iV) 

Vd gd(z) 

{.N-No) \rr. ai + l) 

Vd ^(ir 


,T>Tc 


-T<Tr 


which is in accordance with Ziff. 


(37) 


Now, comparing Eq. (35) and (27) we get, 

P'V^ = - (38) 

X 

For d-dimensional free Bose gas one can obtain from previous equation 

PVd = (39) 

a 

This is an important and familiar relation, PV — when d = 3[6lia[7l[8]. This actually shows equation (38) is a very 
significant relation for the Bose gas irrespective whether they are trapped or free. And in case of trapped bosons effective 
volume and effective pressure plays the same role as volume and pressure in current textbooks and literatures. 


3.6 Heat Capacity at Constant Volume 

Heat capacity at constant volume Cy below and above Ty 


Cv = T{-^)n,v = 


Nk[x{x+l)^gx+iiz)-X^^^] ,T>Ty 
Nkx{x + i)^ax + i) .T<T, 


(40) 
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From this we can investigate whether the specific heat will show a jump or not. From Eq. (40) we can obtain the 
difference between the heat capacities at constant volume, at Tc as 


c 


-Cv 



Nkx^ 


gx(l) 

gx-i(l) 


(41) 


The above equation shows that, the jump of heat capacity depends on y. Cy will be discontinuous at T = if. 



E' 


> 2 


is satisfied. And, Cy will be continuous at T = Tc if, y satisfies 


(42) 


1 < y < 2 (43) 

So, in case of free massive boson (choosing I = 2 and all rn —>■ oo), will itself be discontinuous for d > 4, in agreement 
with Ziff [5] . And in the high temperature limit of Cy approches its classical value as it becomes yAfc for trapped system 
and ^Nk for free system, which is ^Nk, when d = 3. 


3.7 Isothermal Compressibilty 

The Isothermal Compressibilty of trapped Bose gas can be obtained, 

dP' 

QY )n,t 

Now using simple chain rule of partial derivative the above equation becomes 


dP' dP' 

= ( 


dz 


f)z 

)nA-^)n,t 


Using equation (15) and (35), we get 


Which concludes. 


^NkP’ g^[z) 


Avj’ = 


v gy-i(^) 

NkT g^{z) 

v C(x-i) 

NkT C(X) 


,T > Tc 
,T< Tc 


(44) 


(45) 


(46) 


(47) 


which reproduces the same result for isothermal compressibilty of free massive Bose gas at d = 3 [5]. And in the high 
temperature limit kt takes the classical value for free system, which is -p. 


3.8 Heat Capacity at Constant Pressure 

Now, heat capacity at constant pressure Cp, 

Cp = T{^)n,p 

\Nk[{x + lfgl+i{z)g^-i{z){^Y-xix+^)9x+A)^] ,T > 

\ Nklix + imx + l)C(x - 1) W(i4)^ - x(x + l)C(x + l)i4] , T < Tc 
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In case of free massive Bose gas, the above equation reduces to, 


Cp = T{ — )n,p 

r + + ,T>r, 

\ Nk[{i + inli + i)c(i -- f(i + i)c(f + 1 )^] ,t< r. 


(49) 


For d = 3, it coincides exactly with Ref. [SIH] which diverges as j^T' '^ti^n the temperature approaches Tc from above. 
Again in the high temperature limit Cp becomes (y + l)fVfc for trapped system and + l)7Vfc for free system, which is 
|fVfc, when d = 3. So, in the high temperature limit Cp approches its classical value. 


Now, the ratio, 7 = (^) for T > Tc is given by, 




7 = 


9x0 


x{x + 1)i)^t5x+i W - X^^^) 


(50) 


Now, for T >> Tc, the above equation becomes, 

^ (x + 1)^ - x(x + 1) = 1 i 
x(x + 1) - X^ X 

In case of free system, choosing all Ui —> 00 , we get from the above equation, 

, I 

^ = ' + 5 

So, 7 equals |, when d = 3 and I = 2, thus obtaining the classical value at high temperature limit. 


(51) 


(52) 


3.9 Clapaeyron Equation 

In any hrst order phase transition pressure is governed by Clapaeyron equation, in the transition line [3 IS]. Like the 
BEC for free Bose gas at d = 3, BEC for trapped Bose gas will also be a first order phase transition if they obey the 
Clapaeyron equation. The Clapaeyron equation derived from Maxwell’s relations is 


dP As I 
dT Av T Av 


(53) 


where, I, As and Av are the latent heat, change in entropy and change in volume respectively. The effective pressure in 
phase transition line is. 


MT) = 


(54) 


Now, differentiating with respect to T leads. 


dT 


77^[(X+l)fcr 

TVg 


gx+i(i) 

flx(l) 


(55) 


Now when two phases coexist the non condensed phase has specific volume Vg whether the condensed phase has specific 
volume has specific volume 0, which concludes Av = Vg. Which concludes the above equation to be Clapaeyron equation 
where the latent heat of transition per particle in case of trapped boson is 


L = 


gx+i(i) 

gx(i) 


(x + i)fcT 


(56) 


Therefore, BEC for trapped Bosons are also a first order phase transition in arbitrary dimensions. For free massive boson 
the latent heat per partcle becomes, 

agrees exactly with Ref. [7]. 
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4 Discussion 


Properties of ideal Bose gas in the presence of an external generic power law potential (not essentially symmetric) are 
discussed in this section. The study done by density of states approach in d dimensional space allows us to investigate 
the two phase system that is both in the condensed and non condensed phase across the phase transition point. The 
most general conditions for BEC as well as the continuity conditions of specific heat for ideal Bose gas trapped in generic 
power law potential are acquired with any kinematics characteristic in d dimensional space. 

At first, the density of states and grand potential of Bose gas under generic power law potential in d dimension has 
been calculated from which we can easily derive the expression for number of particle, which gives us the first insight 
about the condensed phase. Using number of particle equation, the general criterion for BEC has been obtained from 
which, one can achieve the same conclusion for symmetric potential [17], harmonic potential [18] as well as more simplified 
system of free massive boson [S] [TJIH] . The fraction of condensate as well as the critical temperature was achieved in this 
process, predicting the exact equation for free system]^. Now turning our attention in thermodynamic quantities such 
as, internal energy E, entropy S', free energy A all of which are evaluated from the grand potential. The obtained ex¬ 
pressions coincides with available literature when the potential is symmetric (ni = n 2 = ■■ = rii = .. = Ud)[T7], harmonic 
(ni = n 2 = ■■ = rii = .. = Ud = 2) [18] or the system is free (all rii —> oo)[9]. Now, from the pressure equation we 
have derived the equation of state which determines a relation among the thermodynamic generalized coordinate y and 
the thermodynamic generalized force Y, and temperature T of the system. Point to note in case of trapped boson we 
have seen, pressure P is replaced by effective pressure P' and volume V is replaced by effective volume V. It is seen 
that the high temperature limit of Bose gas is just classical system obeying Maxwell-Boltzmann distribution. Therefore, 
the idea of effective volume and pressure is not only valid for trapped Bose gas but also for trapped classical system as 
well trapped Fermi gas. For free system, both effective pressure and effective volume reduces to pressure and volume. 
So we can say effective volume and pressure to be more general notion for trapped system, which enables us to treat 
the trapped atomic gases as free one. The other derived quatities for trapped system such as specific heat at constant 
volume Cu, specific heat at constant pressure Cp, isothermal compressibilty n, all can reproduce their familiar form found 
in literature [6] |7j | 8 j |9] in case of d dimensional system [ 8 ] [9] as well as in d = 30I2]. 

Now, lets turn our focus towards general criterion of BEC (Eq. 22) and jump of C,, at T = Tc (Eq. 42). In gen¬ 
eral it is seen that there is no first order pi?] IS] as well as no second order [5^ phase transition for d ^ 2 for free system. 
Eq. (25) derived from the general criterion Eq. (22) also suggests the same for the free ideal Bose gas. Now turning our 
attention to trapped system, consider a one dimensional system with ni = 2 and I = 2. Eq. (22) then dictates, no BEC 
for such system as y < 1. In case of two dimensional Bose gas (d = 2), with ni = n 2 = 2 (i.e. harmonic potential) we 
can find it fulfills Eq. (22), so BEC exists in such system. That means, we have seen the case when BEC can exist for 
d ^ 2, with appropiate choice of power law exponent, which was not the case for ideal free Bose gas. In case of jump of 
Cy in free Bose system, it is seen that, there is a jump of C„ at T = Tc, when d > m- For instance, taking d = 3, Z = 2 
and Hi —> oo we see that there is no discontinuity at d = 3 for free system as y < 2. But again one can obtain jump 
in Cy a,t T = Tc, for trapped system when d ^ 4 with appropiate choice of rii. For example, taking d = 3 and I = 2, 
with 111 = 712 = = 2 one can obtain from Eq. (42) that, y > 2. So, in this case we can obtain a jump at d = 3 in 
C„ at T = Tc for trapped system, although not a property for free Bose gas. Neveretheless, Eqs. (40) and (41) shows 
that although Cj.+ is not a monotonical function of y, both Cj,- and ACt^ increase monotonically as the parameter y 
increases. This indicates that, in general, in the vicinity of Tc, larger the value of y, larger the energy needed to raise 
the Bose system to higher temperature state, partticularly when the system is in condensed phase. On the other hand, 
when the temperature is very low, C„(t<Tc) decreases as y increases and the larger the value of y, more quickly C„(t<Tc) 
tends to zero. Also larger value of y indicates larger fraction of bosons to be in the ground state, such that less number 
of bosons excited when temperature raises. 

The Clapaeyron equation applies to any first order phase transition. From the expression of effective pressure it is 
shown that trapped bosons in arbitrary dimension does satisfies the Clapaeyron equation, from which we can say BEC 
is a first order phase transition with or without the existense of external potential. The latent heat of transition per 
particle in case of trapped system is also obtained from this equation, which reduces to exact expression of laten heat in 
case of free massive bosons at d = m- 


5 Conclusion 

From complete thermodynamics of ideal Bose gas trapped in generic power law potential we have derived the general 
criterion for BEC. Also, the calculated physical quantities reduce to expressions available in the literature, with appropiate 
choice of power law exponent. In this manuscript, we have restricted our discussion in case of ideal system. It will be 
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vey interesting to see the effect of interaction on general criterion of BEC. 
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